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SUMMARY 

The birefringence of crystalline substances In general, and of sapphire 
In particular. Is described. A test Is described whose purpose Is to determine 
the direction of the optic axis of a cyllndrlcally machined single crystal of 
sapphire. This test was performed on the NASA Lewis sapphire cylinder and It 
was found that the optic axis made an angle of 18° with the axis of symmetry 
of the cylinder. 


INTRODUCTION 

Amorphous materials such as plexiglass or fused silica and crystals with 
a simple cubic unit cell such as sodium chloride are optically Isotropic. For 
light Incident upon such materials from any direction, the Interaction between 
the Incident light and the material Is described by a single scalar dielectric 
constant K and a single Index of refraction. For all other types of crystal- 
line materials, Incident light from different directions with respect to the 
unit cell of the crystal Interacts with the material In optically unequivalent 
ways. As a result, this Interaction Is described by a diagonal three-by-three 
dielectric constant matrix which Is referenced with respect to a special set 
of axes on the crystal's unit cell. The nonzero elements of this material are 
unequal. Such nonlsotroplc materials are called blrefrlgent or doubly 
refracting. 

For uniaxial crystals (l.e., trigonal, tetragonal, and hexagonal unit 
cells) of which sapphire Is an example, two of the nonzero elements are equal 
and are denoted by K 0 . The third Is denoted by K e . K 0 Is called the 
ordinary dielectric constant and K e Is called the extraordinary dielectric 
constant. Uniaxial crystals possess a single optic axis and all directions In 
the plane perpendicular to the optic axis are optically equivalent. The 
direction of this single optic axis with respect to the crystal's geometry 

will be denoted by N. For the remainder of this report, any symbol with a 
carat over It denotes a unit vector. 

For biaxial crystals (l.e., trlcllnlc, monocllnlc, and orthorhombic unit 
cells), all three diagonal elements of the dielectric constant matrix are 
unequal. Such materials possess two different optic axes and will not be 
considered further In this communication. 
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For uniaxial crystals, the ordinary and extraordinary Indices of 
refraction, n 0 and n e respectively, are the square root of K 0 and 
K e . For any given material, these Indices are weakly wavelength dependent, 
and for sapphire with Incident light of X = 0.546 ym, It Is found (ref. 1) 

that n 0 = 1.771 and n e = 1.763. Sapphire Is said to be a material of 

low birefringence since the difference between n 0 and n e (0.008 In this 
case) Is much smaller than n 0 or n e Itself. 

The standard method for determining the direction of the optic axis of a 
crystal Is to cleave or cut a sample whose thickness Is less than 1 mm and 
examine It between the crossed polarolds of a polarizing microscope (refs. 2 
to 4). However, In a number of situations, It Is Inconvenient or Impossible 

to remove a thin section from the crystal. In these cases, the optical axis 

direction must be determined by using the entire crystal which Is perhaps many 
Inches In diameter. The procedure described In this report represents an 
adaptation of the polarizing microscope procedure to large samples. Such Is 
the case If one Is to perform laser doppler veloclmetry (LDV) measurements 
Inside a firing engine cylinder. Traditional materials such as plexiglass or 
fused silica cannot withstand the large pressures produced during combustion. 
Thus, fabricating optical access windows from high strength materials such as 
sapphire overcomes these mechanical problems. But at the same time, the 
bl ref rlngency of sapphire produces new problems In making reliable optical 
measurements. In order to overcome these optical problems, the direction of 
the optic axis of the access window must be determined without damaging It or 
without removing any material from It (ref. 5). 
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The Geometry of the Ordinary and Extraordinary Rays 

Light waves are coupled electric, and magnetic, it, fields whose 

surfaces of constant phase are perpendicular to the wavevector lc and whose 
direction of energy propagation, the Poyntlng vector, Is given by 

S = t x H /\t x H| (1 .1) 

The restrictions on the possible modes of propagation of light within a material 
due to Maxwell's equations alone leave a substantial amount of freedom as to 
the details of the propagation. All that Maxwell's equations require Is that 
within a material, the wavevector, the electric field, and the electric dis- 
placement vector ft lie In a plane and that ? and ft are related by (ref. 6) 

2 

TU (E - k (k • E)) ( 1 . 2 ) 

where y Is the magnetic permeability of the material and n Is Its Index of 
refraction. The exact details of the propagation of light within a material 
are provided by the constraints Imposed by the constitutive relation between 

E and D, 

^ = [K] £ vac ? (1.3) 

For amorphous materials, the combination of equations (1.2 and 1.3) yield only 

a single solution, light waves for which t Is transverse to £ and for 
which the direction of the wavevector and the direction of energy flow are 
Identical. The phase velocity of such light waves Is 


V 


P 



(1.4) 


and Is Independent of the direction of propagation within the material. 


The combination of equations (1.2) and (1.3) for uniaxial crystals whose 
dielectric constant matrix Is 

"~K 


[K] = 


0 

K 


0 

0 

K_ 


(1.5) 


allows two modes of propagation for the light within the crystal, the ordinary 
and extraordinary modes denoted by the subscripts o and e. The ordinary 
mode Is simple. It Is Identical to the propagation of light within an 

amorphous material. The directions of £ and § are Identical, £ Is 

oo 

transverse to this direction, and the phase velocity of the light Is 
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( 1 . 6 ) 
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independent of the direction of propagation within the material. 

For the extraordinary mode of propagation, the details of the propagation 
vary as a function of the angle <j» between the wavevector k g and the optic 
axis N. Specifically, the electric field vector Is not transverse, 


E • k & * - Q cos sin 


where 


Q = 



(1.7) 


( 1 . 8 ) 


Further, the phase velocity of the light Is direction dependent (ref. 7). 


V 


p 



(1.9) 


A consequence of equation (1.7) Is that the direction of energy propagation 

S Is not Identical to the direction of the wavevector k . This presents a 
e e 

difficulty In that although the mathematical formulation of the behavior of 
the extraordinary mode Is most simply written In terms of k , all experimental 

A ® 

measurements of energy flow measure S g . 


Both the ordinary and extraordinary light waves traveling within a uniaxial 
crystal have definite states of polarization. By solving equations (1.2) and 

(1.3) simultaneously for the components of 6, we obtain the polarization 
directions of the ordinary and extraordinary rays, 


6 = k x N/|k n x N | 

OO 0 


( 1 . 10 ) 


and 

D e = k e x (k e x N)/|k e x (k e x N) | (1.11) 

Solving for the components of the electric field for the extraordinary ray and 
computing S via equation (1.1) we obtain a relationship between S and k 

V C V 

S e - k e - Q(N • k e )N/|k e - Q(N • k e )N| (1.12) 
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This result Is significant since It allows the determination of the experi- 
mentally measured propagation direction using the mathematically convenient 

wavevector k . 

e 

The actual problem at hand Is the following. Consider a uniaxial crystal 
with a flat surface, the outward normal to the surface being n and the direc- 
tion of the optic axis being N. Consider light of an arbitrary polarization 
Incident upon the surface with the wavevector k^ at an angle of Incidence y^ 

as In figure 1. At the surface, the light wave breaks Into two components, 
the ordinary one with a certain polarization direction, and the extraordinary 
one with another polarization. The problem of Interest Is to determine In 
which directions the light waves travel within the crystal and what Is the 
percentage of Initial light energy In each polarization mode. 

The calculation which determines this Information proceeds as follows. 

The continuity of the fields at the crystal surface requires that (ref. 8) 


k k 

crystal _ _1 

U P " C 


a n 


(1.13) 


For the ordinary ray, equation (1.13) determines that the ordinary wavevector 

k has the angle of refraction y as In figure 1, which obeys the Snell's 
o o 

law relation 


sin = n 0 sin y 0 


(1.14) 


and that k. , n, and k = S all lie In the same plane. For the extra- 
1 o o 

ordinary ray, equation (1.13) becomes 


,- 1/2 


2 2 

, /sin $ cos a» . 

s1n Y 1 = ( — 2 + 2 1 s1n 


(1.15) 


where y g Is the angle of refraction of the wavevector k g In figure 2 and 

where k^, n, and k g all lie In the same plane. However, since ¥ Is the 

angle between k g and the optic axis N, It Is a function of y , and as a 

result, equation (1.15) represents a fourth degree polynomial In sin y e . 

The refraction of the extraordinary ray does not obey Snell's law. 

For materials of large birefringence, there Is no choice but to solve 
this fourth order equation numerically. But for low birefringence materials 
such as sapphire where 
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Q = 0.0090 


(1.16) 


an accurate approximation to the solution of the equation may be obtained as 
follows. Let the angle r In the plane of k^ , n, and k g be defined as In 
figure 2 by the Snell's law relation 

sin yi = n e sin r (1.17) 

Then the refraction angle Ye ma y be written In terms of a power series In Q, 


Q sin y. cos 2 e ? 

zn os" 
e 

where © Is the angle between the ray which makes the angle r with the 

surface normal and N. Then once k Is determined from equation (1.18), the 

e 

propagation direction for the extraordinary ray Is given by equation (1.12). 

For normally Incident light, equations (1.14), (1.17), and (1.18) give 

k, = k = k and the polarization directions of the ordinary and extraordinary 
1 o e 

rays of equations (1.10) and (1.11) reduce to 

D q = k^ x N/|k 1 x N | (1.19) 

and 

D e = k 1 x D q / | k^ x D o l (1.20) 

This geometry Is Indicated In figure 3. The geometry of this figure Is 
especially appropriate to the problem of the cylindrical sapphire crystal. 

Let the axis of symmetry of the cylinder be the 2 axis and let N be at an 
angle e with respect to the Z axis In the XZ plane. Let the laser 
beam be normally Incident upon the cylindrical surface In the XY plane at an 
angle of <j> with respect to the X axis. For this case, the unit vector 

A 

°o ,s 


sin <t> cos eu - cos <t> cos eu - 
J; ' y 

(cos 2 e + sin 2 e sin 2 


sin 4 > sin eu z 

T7z 

4 >) 


( 1 . 21 ) 
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and the unit vector D Is 

e 


D = 
e 


-sin <j> sin eu + sin <t> cos $ sin eu - cos eu 
x v i 


2 2 2 
(cos e + sin e sin <t>) 


1/2 


( 1 - 22 ) 


The Determination of the Direction of the Optic Axis 

Consider light normally Incident upon the cylinder which Is polarized at 
an angle a clockwise from the Z axis as Is shown In figure 4. Then the 
direction of the polarization of the Initial light Is 

D 1 = - sin o sin <|>u x + sin a cos 4»u y + cos au z (2.1) 

When this light Is Incident upon a hollow crystalline cylinder, a certain 
fraction A of the Incident light amplitude propagates In the ordinary mode, 
and a certain fraction B propagates In the extraordinary mode as Is shown In 

figure 5. Continuity of the tangential components of E at the Interface 
requires that 


D 1 = AD o + 


BD_ 


( 2 . 2 ) 


Combining this with equations (1.21), (1.22), and (2.1) we obtain 


and 


sin a. cos e » cos a sin <t> sin e 
(cos 2 6 + sin 2 e sin 2 (fr) 1 ^ 2 


sin a sin <t> sin e - cos a cos e 
(cos 2 e + sin 2 e sin 2 4>) 1/2 


(2.3) 


(2.4) 


The light amplitude leaving the front section of the hollow crystalline 
cylinder has the form 


a 


leaving front section 


AD o 


+ e 



(2.5) 


where the e^* term represents the phase delay of the extraordinary ray 
relative to the ordinary ray. This Is due to the fact that within the 
cylinder, the ordinary and extraordinary rays have different wavelengths and 
travel different paths. Thus, the two rays oscillate through a different 
number of cycles and are out of phase when they leave the cylinder, l.e.. 


$ 


2tra 



- %> s1n 
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( 2 . 6 ) 
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For the geometry of figure 3, e reduces to the angle between k^ and N so 
that 


sin 2 © = 1 - sin 2 e cos 2 $ 


(2.7) 


When the light that has passed through the front sector of the hollow 
cylinder traverses the rear sector of the cylinder, that fraction which was In 
the ordinary mode continues solely In the ordinary mode, and that fraction 
which was In the extraordinary mode continues solely In the extraordinary 
mode. The phase difference between the ordinary and extraordinary rays 
Increases from § to 2§. If this light Is then passed through a Polaroid 


crossed with the direction 6^, the resulting normalized transmitted Intensity 
Is 


I 


final 


= [sin 4 i] 


2 2 2 12 
sin 2a (sin 6 sin e - cos e) - 2s1n ± sin e cos e cos 2a f 


2 2 2 
sin <t> sin e + cos e 


( 2 . 8 ) 


The first term In this Intensity expression Is commonly called the Isochromat 
pattern and the second term Is commonly called the Isogyre pattern (ref. 9). 

If the light passing through a thick section of the crystal Is projected onto 
a screen, the Isochromat pattern will appear as a series of from 2 to 20 
fringes across the field of view. The Isogyre pattern will appear as an over- 
all Illumination level or contrast level modulating the Isochromat pattern. 


With this general geometry, the determination of the optic axis proceeds 
as follows. First the Initial polarization 6^ Is set to be vertical, l.e., 
a = 0°. Then the Isogyre contribution to the normalized Intensity becomes 


I 


I 2 sin <i> sin e cos e 
Isogyre ) cos 2 e + s ^ n 2 e s ^ n 2 ^ 


(2.9) 


If the laser source Is stationary and the cylinder Is rotated, then e remains 
constant and <|> rotates with the cylinder. The total Intensity on the viewing 
screen Is zero, l.e., extinction occurs, when <t> = 0° or $ = 180°. This Is 
when the optic axis Is In the plane of the laser beam and the symmetry axis of 
the cylinder as Is shown In figure 6. The Isogyre Intensity contribution 
reaches Its maximum value at sin <t> = ±cot e, or simply at <t> = 90° and 
<t> = 270° when |e|< 45°. Rotation of the cylinder to the maximum Isogyre 
Intensity position produces the maximum bright to dark contrast of the 
Isochromat Intensity pattern which takes the appearance shown In figure 7. 


In this 4 > = 90° or 270° position, with |e|< 45° and variable polari- 
zation direction a, the Isogyre contribution to the normalized Intensity 
becomes 


Ilsogyre = t s ^ n 2(a + ©)] 2 (2.10) 

This has extinctions when a = - e or when a + 6 = 90°. The first 
alternative occurs when the Initial polarization 6^ points In the direction 
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of the optic axis N and the second alternative occurs when the final Polaroid 

which Is crossed with 6^ points In the direction of the optic axis. To 

distinguish between these two alternatives one returns to the configuration 
for the Isochromat Intensity of figure 7, and rotates the cylinder a few 
degrees. This will cause the Isochromat fringe pattern to migrate as In 
figure 8. The direction of this migration Is the direction of the optic axis. 


Experimental Procedure 

The following procedure Is used to align the optical system In order to 
prepare for the determination of the optic axis. 

1. Set the laser In the single line mode and project It onto the viewing 
screen. Mark the laser spot on the screen. This spot will serve as a refer- 
ence location for all further alignments. 

2. Set the cylinder on a rotating table In the laser beam. To ensure 
that the laser beam passes through the center of the cylinder, make sure that 
the light reflected from the front surface of the cylinder appears on the mouth 
of the laser. Also make sure that the laser light transmitted through the 
cylinder Is centered on the marked spot on the viewing screen. 

3. Set the microscope objective to be used as a beam expander Immediately 
In front of the mouth of the laser. The laser light on the viewing screen 

will now be very diffuse and the Illumination region will have a large diameter. 
Center It as best as Is possible on the marked spot. 

4. Set the beam stop with a roughly 1/8 In. diameter hole between the 
microscope objective and the cylinder. The brightest center portion of the 
light from the microscope objective should pass through the hole In the beam 
stop. The resulting light pattern on the viewing screen should be centered on 
the marked spot. This step In the procedure serves to filter the laser beam. 

5. Set the polarization rotator between the microscope objective and the 
beam stop near the focal point of the objective lens. The light pattern on 
the viewing screen should be centered on the marked spot. Set the polarization 
to be vertical. This vertical alignment Is Important because, as will be seen 
In the next section, the location of the optic axis Is very sensitive to mis- 
alignments In the Initial polarization. At this point note also that for many 
polarization rotators every 1° rotation of the polarization rotator dial pro- 
duced a 2° rotation of the polarization of the laser light. It Is also Impor- 
tant to note that If a series of mirrors are used to translate the laser beam, 
the polarization should be aligned after the last mirror. This Is because 
when polarized light diagonally reflects off of a mirror, the direction of the 
polarization of the light Is changed. 

6. Set a Polaroid between the cylinder and the viewing screen. The Pola- 
roid should be crossed with the polarization rotator. This Is accomplished by 
rotating the polarold until the fringe pattern on the viewing screen has the 
poorest bright to dark contrast. It Is Important not to place the polarold In 
an unexpanded laser beam. If this Is done, some types of polarolds burn out 
In about 2 sec. 
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This completes the alignment of the optical system, and the geometry of 
steps 1 to 6 Is shown In figure 9. At this point, the direction of the optic 
axis Is determined In the following way. 

7. Observe the Interference fringe pattern and its bright to dark contrast 
on the viewing screen. Rotate the cylinder until extinction occurs. The 
extinction lasts for a rotation range of about 5°. Choose the center of this 
range as the $ = 0° position. The plane of the laser beam and the axis of 
symmetry of the cylinder now contains the optic axis.’ 

8. Rotate the cylinder 90° until the resulting fringe pattern appears as 
In figure 7. Notice the way the fringes migrate under small rotations of the 
cylinder as In figure 8. Now In the 90° position, rotate the polarization 
rotator and the crossed Polaroid until extinction occurs. The polarization of 
the Incident laser light Is now parallel or perpendicular to the optic axis. 

The direction of the fringe migration of figure 8 enables one to choose the 
direction of the optic axis from these two alternatives. 

When this procedure was performed on the NASA Lewis sapphire cylinder, It 
was found that the optic axis made an angle of roughly 17° to 19° to the ver- 
tical as Is shown In figure 10. This cylinder has an outer radius of 2 3/8 In. 
and a wall thickness of 5/8 In. 


Polarization and Stability Questions 

Referring to the geometry of figure 3, If the optic axis were along the 

symmetry axis of the cylinder, 0 g would always be In the Z direction and 6 q 

would always be In the XY plane normal to the laser beam direction for all 
angles <t>. Thus, vertically polarized Incident light would propagate solely 
In the extraordinary mode and horizontally polarized Incident light would propa- 
gate solely In the ordinary mode. As a result. In the crossed Polaroid experi- 
ment of Section 4, extinction would occur for all rotations of the cylinder <t> 
for vertically or horizontally polarized light. Similarly, If the optic axis 

of the cylinder were at e = 90° along the X axis as In figure 3, D g would 

always be In the XY plane normal to the laser beam direction and D q would 

always be along the Z axis for all angles $. Thus, vertically polarized 
Incident light would propagate solely In the ordinary mode and horizontally 
polarized Incident light would propagate solely in the extraordinary mode. 

Again for the crossed Polaroid experiment, extinction would occur for all 
rotations of the cylinder <t> for vertically or horizontally polarized light. 

For the situation at hand when 6 Is neither 0° nor 90° the polarization 
directions are In general not as simple as for the two previous cases. The 
only case where simplicity arises Is when the cylinder Is rotated to the <|> = 0° 
position as In step 7 of the previous section. In this case, the plane of the 
laser beam and the symmetry axis of the cylinder contains the optic axis as In 

figure 6. For this geometry, D g Is along the Z axis and D q Is along the Y 

axis. Thus vertically polarized Incident light propagates solely In the extra- 
ordinary mode and horizontally polarized Incident light propagates solely In 
the ordinary mode. For this geometry, although birefringence Is possible. 
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with vertically or horizontally polarized light It In fact does not occur. It 
Is expected that this geometry will be of central Importance to LOV applica- 
tions. 

In the previous section, Instabilities due to Initial polarization 
misalignments were alluded to. These are calculated In the following way. 

For the angles e, <t>, and a In figure 3 being arbitrary, the Isogyre inten- 
sity Is given by equation (2.8). Extinction occurs when that expression Is 
equated to zero, l.e., when 


sin <t> = 


tan a 
tan e 


(3.1) 


As seen In the previous section, for vertical polarization, a = 0° and 
thus the <t> * 0° plane which contains the optic axis Is determined. For 
small misalignments a and with 6 = 18° for the NASA Lewis cylinder, for 
every degree that a Is misaligned, the value of $ where extinction occurs 
Is roughly 3° away from the <t> = 0° position. Thus any misalignment In the 
Initial polarization Is magnified by a factor of about three In the determina- 
tion of the direction of the optic axis. 


CONCLUSIONS 

For small, thin samples of a blrefrlngent material, the polarizing micro- 
scope Is a useful tool for determining the direction of the optic axis. For 
large samples, a variant of the polarizing microscope procedure has been des- 
cribed which allows one to Identify the optic axis direction to within ±1°. 
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APPENDIX 


An Interesting consequence of equation (1.12) Is that although the ordin- 
ary ray always propagates In the plane of Incidence within the crystal, the 
extraordinary ray almost never does. This Is especially noticeable for the 
case of normal Incidence. For normal Incidence, the ordinary ray continues to 
propagate straight ahead undeflected while the extraordinary ray Is deflected 
at the surface of the crystal towards, or directly away from the direction of 
the optic axis. As long as the ordinary and extraordinary rays propagate 
within the crystal, the distance between them Increases. However, once they 
leave the crystal through another flat surface parallel to the entrance surface, 

they propagate parallel to each other In the direction of with the separ- 

ation f. For normal Incidence, the distance f Is calculated as follows. 

Let 0 be the angle between and N. Let % be the angle between S g 

and k^. Then equation (1.11) gives 


cos \ 


1 - Q cos 2 0 

[1 + (Q 2 - 2Q) cos 2 ©] 1/2 


(A.l ) 


For small Q, the denominator may be expanded to give 


sin 


1 a 

2 ~ 4 


sin 20 


(A. 2) 


Finally, the separation f Is 


f 


a tan % ~ 


9i 

2 


sin 20 


(A. 3) 


For normal Incidence through a 1 -1 n . thickness of sapphire 


f < (1 In.) (25 400 ym/1n.) 


0.0090 


= 114.2 ym 


(A. 5) 
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Figure 1. - The geometry for the refracted ordi- 
nary light ray produced at the interface of a 
uniaxial crystal 



Figure 2. - The geometry for the refracted extra- 
ordinary light ray produced at the interface of 
a uniaxial crystal. 


Figure 3. - The geometry for light normally incident upon a 
hollow cylindrical crystal and the ordinary and extra- 
ordinary polarization directions within the crystal. 



Figure 4. - The initial polarization of 
the normally incident light ray. 



Front sector of 
the hollow 
cylindrical 
crystal 

Figure 5. - The accumulation of phase delay of the extraordinary ray 
relative to the ordinary ray. 



Figure 6. - The directions of the ordinary and extraordinary ray polarizations 
when the optic axis is in the plane of incidence. 




Figure 7. - The isochromat inten- 
sity pattern when the optic axis 
is perpendicular to the plane of 
incidence. 



Figure 8. - The migration of the isochromat inten- 
sity pattern. 
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Figure 9. - The optical system schematic. 
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Figure 10. - The direction of the optic axis for the 
sapphire cylinder used by the NASA Intermittent 
Combustion Engine Branch. 
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